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^ ' Abstract 

. We consider the stability of the two branches of non-extremal enhangon solu- 

tions. We argue that one would expect a transition between the two branches at 
some value of the non-extremality, which should manifest itself in some instabil- 
ity. We study small perturbations of these solutions, constructing a sufficiently 
OO \ general ansatz for linearised perturbations of the non-extremal solutions, and 

' show that the linearised equations are consistent. We show that the simplest 

Q ■ kind of perturbation does not lead to any instability. We reduce the problem 

CN , of studying the more general spherically symmetric perturbation to solving a 

^ ' set of three coupled second-order differential equations. 

1 Introduction 

^ ! A key issue in string theory is the role and physical interpretation of singularities in 

\ supergravity solutions. Some singular solutions, such as negative mass Schwarzschild, 

\-{ ' are genuinely unphysical [0, and are simply excluded from consideration; no corre- 

■ — ' sponding source exists. String theory provides resolutions of many other singularities 

through various mechanisms. Recently, new singularity resolution mechanisms have 
played an important part in the understanding of field theories with partially broken 
supersymmetry in the AdS/CFT correspondence [0, ^ ||, ||, ^ |^. A simple example 
of this new class of mechanisms is the resolution of the repulson singularity of Q 
by the enhangon mechanism |]10 . 



Generally, one of the simplest questions to consider from the bulk spacetime side of 
the AdS / CFT correspondence is the finite-temperature behaviour of the theory. One 
would expect that the theories with reduced supersymmetry should have interesting 
phase structures. At high temperatures, one would expect to find that the partition 
function is dominated by a black hole solution, and there may be some symmetry- 
breaking phase transitions as the temperature decreases. Attempts to investigate 
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these issues by studying black hole solutions on the AdS side were made in [11, 12 



13| , |TJ, |T6|. Considerable progress was made on obtaining suitable black hole 
solutions. However, because of the complexity of the setup, no exact closed-form 
solutions are available. 

In this paper, we will begin an investigation of the phase structure associated 



with non-extremal black hole versions of the enhangon solution of [10], using the 
simple explicit solutions generalising the enhangon found in [|Ty, |T^. We will focus 
on studying whether the solutions have classical instabilities which could provide the 
mechanism for transitions between them. We analyse the linearised perturbation 
equations around the non-extremal enhangon background, generalising the analysis 



of [|T8[ in the extreme case. Although the enhangon is somewhat different from the 
asymptotically AdS cases, the underlying physics should be similar. It would be 
interesting to extend our work to consider the stability of the non-extremal fractional 



brane solutions of [jT9[, which are more closely related to asymptotically AdS cases. 

In section |^, we review the extremal and non-extremal enhangon solutions. There 
are two branches of non-extremal solutions, arising from an ambiguity of a choice of 
sign in the solution of the supergravity equations. One branch joins on to the ex- 
tremal enhangon solution studied previously, and always has a shell of branes outside 
the horizon. (The proportion of the energy carried by the shell and by the black hole 
inside the shell in this solution was not determined at this supergravity level; a better 
understanding of the internal dynamics of the shell is required to obtain a unique 
solution for given asymptotic charges.) The other branch appears at a finite value of 
the non-extremality parameter. Above this critical value of the non-extremality, both 
types of solution are possible. At large energies, the effects of the D-brane charges 
should be negligible, so the solution with a horizon, which for large mass is approxi- 
mately the usual uncharged black hole solution, has the right physical behaviour. On 
the other hand, if we begin slowly adding energy to the extremal enhangon, we will 
obtain a solution on the branch with a shell. We would expect that there is some 
non-trivial transition between these two families of solutions as we vary the energy.[| 

We are going to focus on the hnearised stability analysis, but we will begin by 
discussing the thermodynamic aspects.0 In section ^ we will compare the entropies of 
the two solutions, and see that the horizon branch has larger entropy at large mass, 
as we would expect. We can also calculate the specific heat for the horizon branch; 
for the branch with a shell, the ambiguity in the division of energy between the shell 
and the black hole prevent us from obtaining a well-defined answer for the specific 
heat. 



^ Since the enhangon is hke a monopole solution, we expect the physics to be similar to that of 
the Einstein- Yang Mills Higgs system (see and references therein). For any given value of the 
asymptotic charges, only one of the two solutions should be stable. However, to see this physics, it 
may be necessary to include the effects of the non-Abelian gauge fields, as in which we do not 
do. 

^One interesting suggestion in |ll| was that in some cases, black hole solutions should exist only 
for temperatures greater than a critical value. We will see that for the non-extremal enhangon, 
solutions with a regular event horizon exist for only for values of the non-extremality parameter 
greater than a critical value — that is, for sufficiently large energies. There also appears to be a 
maximum temperature for these solutions, but no minimum. 
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Our main focus is to look for dynamical instabilities that could take us from one 
branch to the other. We particularly want to see whether there is an instability at 
some value of the energy which could take us from the branch with a shell to the 
horizon branch, which we think should be the physical solution at large energies. 

In section |, we set up an appropriate ansatz for the perturbations. We consider 
only perturbations which are spherically symmetric in the transverse space and trans- 
lationally invariant along the branes, as we are looking for a transition between two 
solutions which preserve these symmetries. We consider the most general ansatz con- 
sistent with the assumed symmetries. This ansatz is slightly more general than the 



ansatz for perturbations of the extreme enhangon considered in [|T8|; we find that our 
more general ansatz is necessary to obtain non-trivial solutions of the full set of field 
equations. We use the remaining diffeomorphism freedom to reduce the linearised 
equations of motion to four second-order equations for four functions characterising 
the perturbation. One of these equations is decoupled from the others. 

In section we consider the stability to this decoupled mode. This equation is 
in fact identical to the free scalar wave equation in this background. Since the mode 
is not coupled to the shell, it satisfies simple continuous matching conditions there. 
We reduce the equation to a one-dimensional bound state problem, and find that the 
potential is negative in a region near the shell, so one might expect that there could be 
bound states (and hence an instability). Nevertheless, we present a general argument 
that there can never be an instability associated with this mode. The idea is that 
since the equation is just the free wave equation, a constant function is a solution. 
This implies that the bound state problem has a zero-energy solution with no nodes, 
and as a consequence, there can be no bound states. 

In this paper, we will not consider the solution of the other three coupled equa- 
tions. The boundary conditions at the shell will be more complicated for these modes, 
and we will need to solve the equations numerically to determine if there is any in- 



stability. This analysis will be continued in a forthcoming paper |22 



2 The enhancon solutions 

The original repulson geometry ^ is constructed by wrapping D(p-|- 4)-branes 
on a K3 manifold of volume V. We will also consider including M Dp-branes parallel 
to the noncompact directions of the D(p -|- 4)-branes. This leaves an unwrapped 
l)-dimensional worldvolume in the six non-compact dimensions. There are 5 — p 
non-compact spatial dimensions transverse to the branes. We will consider the case 
p = 2, so we have coordinates r, 6, in the transverse directions, and x'^, /i = 0, 1, 2 
in the noncompact directions along the branes. The Einstein frame metric and fields 
are 

ds' = Z^'/'Z,'/\,dx^dx'' + zTzy\dr^ + r^dn) + V^'^zT Z.^'^'dsl^ , 
e^* = glZ,'/'Z,^"\ 
C(3) = {Z2gs)~^dx^ A dx^ A dx^ , 

C(7) = {Zf^gsY^dx^ ^dx^ ^dx^ Ekz , (1) 
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where the harmonic functions are 



1 + ^, ^2 = 1-^, (2) 



the parameters are related by 



r6 = ^^, -2 = -^re (^1 - - j , (3) 
and dVL denotes the metric on the unit two-sphere. The running K3 volume is 



V = V^ = (27rvQ;') at the enhangon radius, 



e 



-^6 • (5) 



V-V, "V 2NJ 



The repulson singularity would occur aX r = r2 < r^. 

The enhangon mechanism discovered in resolves this repulson singularity. 
The essence of the mechanism is that the singularity can never be formed. If one 
tries to assemble the repulson from well-separated branes, the constituent branes 
will stop behaving as pointlike objects and smear out into extended solitons at a 
certain distance from the would-be singularity; the sphere at this radius is called the 
enhangon sphere. This effect is due to the appearance of additional light degrees of 
freedom, enhancing the gauge symmetry from U{1) to SU{2), when the K3 volume 
is V = = {271 y/a')'^. The metric outside the enhangon sphere is still the repulson 
geometry, but the sources are distributed over the sphere, leaving flat space inside 
and removing the singularity. This is the enhangon geometry. 

Although this singularity resolution depends on stringy physics, namely the ap- 
pearance of additional light degrees of freedom which are not contained in the original 
supergravity description, it was found in that the appearance of a shell at the en- 



hangon radius can be understood from a purely supergravity argument. If we imagine 
distributing the sources on a spherically symmetric shell, so that the exterior geome- 
try is the repulson, while the spacetime inside the shell is fiat, then the energy density 
of the shell will be positive only if the shell lies outside the enhangon radius. Thus, 
the enhangon radius provides a minimum position for the shell. 

Thus, the above geometry does not always apply for all r. For M > N there is 
no repulson singularity, and we can assemble sources to form the geometry in 
For M < N however, some of the D6-branes cease to be pointlike before we reach the 
singularity at r = r2, and will form an enhangon shell. This geometry then applies 
only outside the shell. 



^For M > 2N, there is no enhangon, and we can form the above geometry by bringing the branes 
in individually from infinity. For N < M < 2N, D6-branes on their own smear out at r = ro. We 
can still form the geometry (0) if we first form D2-D6 bound states, which can be brought to the 
origin. 
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We will mainly be interested in M < A^. We assume that all M D2-branes 
coincide at the origin, along with A^' D6-branes, where A^' < M. The remaining 
N — N' D6-branes lie on an enhangon shell. The geometry inside the shell is 

gl'^ ds' = H^'^'H^'/\,dx^dx'' + Hl'^'Hl'^dr^ + r^dn) + V"^hT H^'^'dsl^ (6) 

and the non-trivial fields are 

C(3) = {9sii'i)~^dx^ A dx^ A dx^ , 

C(7) = {gsE^Y^dx^ Ndx^ Ndx^ NV e^^'i , (7) 

where 

H2 = l- - ^, r = re-^— ^ , 8 

Te r V IS 

H, = l + ^ '- + ^, r=r6— = ^^- . 9 

Te r I\ 2 

The constant terms in the harmonic functions are chosen to ensure continuity of the 
solution at the shell. 

The supergravity argument can be extended to non-extremal generalisations of 
the enhangon solution, which are difficult to study from the string theory point of 
view. A non-extremal solution was first written down in |T^. In it was found 
that there are two branches of non-extremal solutions, arising from an ambiguity of a 
choice of sign in the solution of the supergravity equations for the usual ansatz. The 
non-extremal generalisation of the exterior geometry is 

gl/^ds' = Z^'/'Z,'^\-Kdf + dxl + dxl) + Z'^/'zl^\K-'dr^ + r^dQl) 

+V'/^ZrZ^'/'dsl, , (10) 

the dilaton and R-R fields are 

^2$ _ 27I/27-3/2 

C(3) = {gs<^2Z2)^^dt A dx^ A dx"^ , 

C(7) = (gsaeZQy^dt A dx^ A dx'^ A V eK3 , (H) 
and the various harmonic functions are given by 

K = 1-!°, 
r 

Z2 = l + - ^6 = 1 + ^. (12) 



r 



Here 



h = -i- + \ri+i^] , (13) 
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and ag = ^&/^&- There are two choices for f2 consistent with the equations of motion: 




^2 = -^±tri+h^ , (14) 



and a2 = T2/r2- Here, r2 and rg are still given by (^. We have changed our conven- 
tions for Z2 to facilitate comparison of formulae involving Z2 and Zg, so the repulson 
singularity, if there is one, is at r = —r2- 

There are two branches of solutions. For the upper sign in (|14|), f2 > 0, so there 
is no repulson singularity, and the solution has a regular horizon at r = tq. For 
the lower choice of sign, however, the repulson singularity always lies outside the 
would-be horizon, \f2\ > tq, and the geometry will be corrected by an enhangon 
shell. We therefore refer to the former as the 'horizon branch' and the latter as the 
'shell branch'. It is interesting that the appearance of a repulson singularity in the 
non-extremal solutions is not connected to whether M > N, but rather to a discrete 
choice. For M > N, the extremal solution is the same as the solution at ro = on 
the horizon branch, and we regard the horizon branch as the only physical choice. 
For M < N, on the other hand, the extremal solution is the solution at = on the 
shell branch, so we need to consider both branches of solutions. We will henceforth 
focus on the case where M < N. 

The shell branch exterior solution is cut off by an enhangon shell at 

_ Krg - Vh 

As in the extremal case, this shell will contain N — N' D6-branes, while the interior 
solution with M D2-branes and N' < M D6-branes is 



H-'^H-'" (-^Ldt' + dxl + dxl] + HTHl'\L-Hr^ + r^Q) 

+V^I^HrH,'"dsl, , (16) 



with accompanying fields 



,2$ _ 2rrl/2rr-3/2 



C(3) = -jj-^ {gsa'2H2)~^dt A dx^ A dx^ 




1/2 

,2 



1/2 

C(7) = [ "'^'"^ ) {gsa'^H^y^dt Adx^ Adx^ AV Skz , (17) 



where 



1 - 



r 



r 



To — K f'n 
H2 = 1+- - + - 
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+ 



\ 



"2 



re 



V, M-N' 
V N 



Ha 



1 + 



re 



r 



-- + 
2 



r';?+r4 



If 



and a'g = r'^/rQ, 02 = ^2/'"2- Note that we have introduced an independent non- 
extremahty scale rg for the interior solution. Implicitly rg < rg in order that the 
interior black hole actually fits inside the shell. We have taken the horizon branch 
for the interior solution, as A^' < M. 

The shell branch solutions have an additional parameter, Tq, which is not deter- 



mined by the asymptotic charges of the solution. It was argued in ||T^ that this was 
simply a weakness of the supergravity excision procedure, and that a better under- 
standing of the physics of the shell should fix this parameter. We will not attempt 
to resolve this issue in this paper, but will simply consider the stability of the shell 
branch solutions for arbitrary Tq < tq, Tg. 



3 Thermodynamics 

We would like to briefly compare the behaviours of the two branches. The ADM 
energy density for these solutions is 

where G is Newton's constant. For the horizon branch, this gives 



EH, = ^iro + \r4 + rl + J'4 + ri\, (20) 



4G \ V 4 " V 4 

while for the shell branch, 



The difference between the ro = solutions is AE = \r2\/2G. For M < A^, we need 
to add this much energy to the extremal solution before we can get solutions on the 
horizon branch. 

The entropy and temperature on the horizon branch are easily obtained from the 
metric (p!OD, giving 

Sn, = ^ = ^{ro + f2f\r, + r,)y^ (22) 
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'-hb - T—r_ — , _^ M/9/._ — , _^ M /o • (23) 



4vr(ro + f2)V2(ro + r6)V2- 



For the shell branch, we must use the interior solution (|T^), which gives 



Ssb = ^ = ^H,{r',f''H,{r',r\ (24) 




" 4^^ U(re)i^2(r^,)if6(r(,) ' ' ^^^^ 

On the horizon branch, we see that the temperature is a monotonic function of rg, and 
hence the specific heat is always negative.^ For the shell branch, we cannot evaluate 
the specific heat, as we do not know rQ(ro). 

The ambiguity in the interior solution on the shell branch prevents us from com- 
paring the entropies of the two solutions for most values of the parameters. However, 
we can make a comparison at large energies, when tq ^ T2irQ. Then 

2 

Tf\ ITT 

Ehb^^, Sf,, ^ ^ 4nGEl, (26) 
as for an uncharged black hole, while for the shell branch. 

Since t'q < tq and V^/V is a small parameter, we conclude that the entropy is larger 
on the horizon branch at large mass. Thus, at least for large fixed mass, we would 
expect the horizon branch to dominate. 

It would also be interesting to compare the entropies at fixed low temperature 
(so again tq ^ r2,rQ). Unfortunately, this is not so straightforward. On the horizon 
branch, 

^ Shb - (28) 



but on the shell branch. 




AnrQ 



(29) 



so 



Thus, whether Sgb is smaller or larger than Shb in this regime depends on how close 
v'q can be to tq. Surprisingly, if it is sufficiently close, Sgb can be the larger. 

*As a consequence, the conjecture of presumably implies that if the xi,X2 directions are 
non-compact, this solution has a Gregory-Laflamme type instability ||2^ . This is not the instability 
we are interested in considering, as it seems unlikely to mediate a transition between our two families 
of translationally-invariant solutions. 
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Thus, we see that thermodynamic considerations suggest that at least for large 
masses, the horizon branch should be the preferred one. Detailed investigation of 
the thermodynamics is hampered by the fact that we don't know how r'^ varies with 
tq. Black hole thermodynamics depends on studying the static vacuum solutions as 
functions of the parameters, so the presence of an unphysical one-parameter ambiguity 
in our family of solutions is a serious impediment. 



4 Perturbation ansatz 



We now turn to our main objective, the consideration of the stability of these so- 
lutions. We wish to consider the simplest set of linearised perturbations of the en- 
hangon solutions which could provoke a transition between the two branches. We will 
therefore assume that the perturbations preserve many of the symmetries of the back- 
ground (|ToD . Specifically, we assume the spherical symmetry in the 6,(f) directions, 
translational invariance in xi and X2, and the discrete symmetries under xi —xi, 
X2 — >■ —X2, 4> — > —4> are preserved. By a suitable choice of coordinates, the most 
general perturbation consistent with these symmetries can be written as the metric 



dilaton 

and R-R fields 
Here 



-V'i/2 



+Zl'^Zl'\K-^dr^ + r^d^l) 



V^'^Zl"z^"^ds 



+ 



C(3) = C'CS) + <5C'(3), C(7) = C(7) + <5C(7)- 



V-i = (/) + Z2 = ^2(1 + SZ2), Zq = Ze{l + 6Ze), K = K{1 + 6K), 



(31) 

(32) 
(33) 

(34) 



the harmonic functions Z2,Zq,K are as in (|1^), the unperturbed dilaton is as in 
(p!T|), and the R-R potentials are as in (|ll|). The first-order perturbations are all 
functions of (r, t) only, while the background quantities are functions only of r. We 
look for perturbations of the form /(r)e*'^*. 

Our ansatz is slightly more general than the ansatz adopted in the study of per- 
turbations of the extremal enhangon geometry in |T^. We have introduced three 
new perturbation functions, 61^2, Sips, and 6K. As we will see shortly, we can choose 
to set SK = by a gauge transformation. The first-order function 61^3 is the only 
thing that breaks the rotational symmetry between xi and X2- As a consequence, it 
decouples from the other perturbations. We could set it to zero without affecting the 
other modes; instead, we retain it, and study it independently of the others. This 
provides us with a single simple (but non-trivial) perturbation equation, which we 
study in section |^. One might think that 6ip2 would also decouple, as it breaks the 
boost symmetry between t and Xi,X2 which ( pUj ) respects. However, the assumption 
that the perturbations depend on t and not on Xi,X2 also breaks this symmetry, so 
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we will find that 61^2 couples to the time derivatives of the other perturbations, and 
it is not possible to set it to zero. That is, it is necessary to consider the more general 
ansatz containing 6ip2 to satisfy all the field equations, even in the extreme case. 

We will now consider the full set of linearised equations for the perturbations. 
The gauge field equations give 

drSC(^S) = -\{^SZ2 + 50 - 5^l)9,C(3) (35) 

and 

dr5C(j) = -^{A6Ze - 36(f) + 35^i)9,C(7). (36) 

The linear part of the stress tensor only involves drSC(^3^ and dr5C(j), so we can 
substitute (|S5|JH0) directly into the stress tensor. 



There are seven distinct equations coming from the linearised Einstein's equations: 
six different diagonal components, and an off-diagonal \tr\ component. With the 
dilaton equation, this gives us eight equations; but with the gauge field perturbations 
fixed by (^,0), there are only seven undetermined functions in our ansatz. The 
problem seems overdetermined, so it is important to ask whether there will be any 
non-trivial solutions of the full set of equations. We have written down the most 
general perturbation consistent with the symmetries we have assumed, so we expect 
there is sufficient redundancy in the equations to admit non-trivial solutions. 

In fact, we can see directly that there are non-trivial solutions to these equations, 
using a trick from |^5|]. We observe that the ansatz (|31|) does not completely fix the 
gauge, as there are infinitesimal diffeomorphisms which preserve its form. Namely, 

t-.t' = t + e*'^*5t(r), r ^r' = r + e''^^5r{r), (37) 

with 

dr5t = iiu^^dr. (38) 

If we apply this diffeomorphism to the non-extremal enhangon geometry (|T0|), we 
obtain a metric of the form (|3ll) with 

<5^f = (^(f)' -^yr-'^dr6r-'^tuj6t, (39) 

. 4 4 4 
Sip2 = Sr + -drSr + -iuj6t 



3r' ' 3" ' 3 ' 




( Z' l\ 2 2 
= M + — Ur - -dM - -iuj5t, 



Zo 3r ) 3 3 



3K'^ = (^ + -]sr- 2dr6r, 
\ K r I 



Ld 



(j)'5r. 
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Since this particular perturbation comes from a coordinate transformation, it must 
solve the equations of motion. Thus, there are non-trivial solutions of these equations. 
Of course, we are not interested in solutions which are pure gauge, but this serves to 
demonstrate that there is some redundancy in the equations. 

This diffeomorphism contains an arbitrary function; since we are not interested 
in pure gauge perturbations, we should fix this additional gauge symmetry. We can 
do so by setting one of the perturbations to zero. It is convenient to choose 5K = 0. 
There remain diffeomorphisms which will preserve 6K = 0: these have 



and 



St 



luja 



6r 



2(ro + f2)('^o + ^6) 



(40) 



+ 



15rQ 



+ 



K{r) 
3ro(f2 + fe) 



+ 



^ + ^ + r2 + h]r^K{r)+ (41) 



f2r6 In 



1 



Kir) 



iub 



where a and b are arbitrary constants. The perturbations (|39| ) with this St, Sr 
then give a two-parameter family of solutions of the linearised equations with SK = 0. 
We will exploit this remaining coordinate freedom to simplify the equations later. 

Having set SK = 0, the contributions to the Ricci tensor linear in the perturbations 

are 

^{2Si>2 + 9Si^i - 5SZ2 + 3SZe) 

^16 ( Sij'i + -Si)'^ - 8 f 5< + -Sij'^ - 8 f SZ'i + -SZ'^ 
\ r J \ r J \ r , 



SR, 



(42) 



32Z2ZQ 

-8 [szi + '^-SZ'^ + m^'^^ + 4#i 



2 
— ( 

r 

^K' Z'^ Z'^' 

40 h 5— + — 

K Z2 Z,, 



-SZ' 5 



Z' 



\ Z2 

+ {Slp2 - SZ2 



z^ 



SZ' -32 



K^ 
~K 



15^ + 3^ 
^2 



Z| K Z2 Zl KZ,, 



SR 



tt + \k {SRu + SR22) = - (5^2 + 85^1 - 65^2 + 2SZ, 



(43) 



K^ 



32Z2Z, 



/ 9 \ /<'' 

24 ( (5V2 + -'5^2 ) + —(245^2 - 325^; - 245Z' + %SZ'^ 
\ r J K 



K Za 



K {SRn - SR22) = Si^s + 



K^ 



8Z2ZQ 



iSip'^— + Sips 



7/2 



K' Z'^ 

Zl ^ ^1{Y2 



+ -Stp'r^ 

Zl K^Z^ 
Zl K Z. 



(44) 
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SR 



tr 



K \ K Z2 Zq 



^\ K Z2 Z. 



+ SZ. 



K' Z' Z'' 

-4 — + 3— 

K Z2 Z. 



(45) 



■^2-^6 

Ho 

1 



+ 5Z2 + 5Z^ 



+ 



32 



K' 



mip'i - 246z!^ + mz'^ - mip'^— + ^v^i 



16 



K' 



Z' 



5Re 



f^2Z& 

4K 
r^K 



(-5^1 + SZ2 + SZ, 



(47) 



+ 



32 



8S1P'; - 85Z'^ - 85Z'^ + 5^P[ + 8^ + 12^ + 28^^ 



'32 K' Z' Z'" 



and 

SRmn[K3] 



VZ2 



4K 



(-5V^i + 5Z2 - 5Ze) 
2 



(4^ 



32Zfi 



8 + -6ij[^ - 8 [dZ'i + -5Z'^ + 8 [sz'l + -5Z'^ 



K' 



. Zr, . Za 



K' 



Zr, Za 



8— + 12^ -4^]+SZU -8— - 3^ + 



K 

K' Z' 



K 



+SZ' 8— + 9-^-3 



Z'' 



K 



Z2 Zq^ 

K Z2 



712 



K^Z^ 
K Z. 



where we have introduced certain combinations which simphfy the resulting equations, 
■ signifies dt, and ' signifies dj.. 

The hnearized Einstein's equations give seven equations. First, there are the 
simple equations from 5_ and which are respectively 







(49) 



and 



^2^6 ^^^^ ^ g^^^ _ g^^^ ^ 25^6) + 3(5^^' + U^P'^ (1 + ^) + (-4#i - 36Z'e +SZ!,)^ = 0. 



,2 K\ 
r 'K 



(50) 
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There are three more independent second-order equations, 




r K Z, 



6. 




- 2^(4(5^^1 + 25i^2 + 6Z,) + Q{5Z'i + -5Z'^) + 5^l^[ - 12^] (52) 

K"^ r \ K Z2 J 



and 

- 8^{7Si^, + 25^2 - + 5Z^) + 245< + 5^l:[ (— + 40^ - 12^ + 36^] (53) 
ii"' \ r K Z2 Zq J 

, , / 24 ir' Z' Z'\ , , /72 /T' Z'\ 

+3(-5(/) + 55^|Jl - 55Z2 + 35Ze) " 
+Q(^_SS<f>-5,lj, + 5Z2 + 5Ze)(^^-^^^ = 0. 

The remaining Einstein's equations give us two equations which are first-order in dr- 
Integrating the tr equation in t gives 

4#^ + 16#; - 86Z'^ + 86Z'^ - 2Sij2^ + #1 (-8^ + 5^ + ^] (54) 

A \ Z2 ZqJ 

-^^^(4-|-|)+^^«(-4-^|-|)+^H-|-^l) - 

and a suitable combination gives the equation 
16^ (4#i + <5V;2 - 2SZ2 + 25Z,) - 8Sf,^ + - 32^ - 12^ - 28^") (55) 



^ K 


+ #'1 


128 

r 


-32^ 


^2 

- 12^ 
Z2 


4^1 




^ 96 
V ^ 


-16^: 


-12^ 
^2 



-Kf -4 -I -I) -4 -I -I) 
(-^f ^ ^^1) + (*^' - (^"i - ^41 + 4 - 41) 
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Finally, there is the dilaton equation 



K'Z'A 





+3(350 + 6^, -6Z,-6Z,){^-—^]+ (50 - 55^1 + 55^2 - 35^6; j ^ ^ ^ 

These equations are coupled in a complicated fashion, but we see that as mentioned 
earlier, there is one simple equation, (|49|). In fact, this is the free wave equation. We 
will discuss the analysis of this decoupled mode in detail in section ^ 

To simplify the other equations, we exploit the remaining two-parameter family of 
diffeomorphisms (|40| , |41|) . These can be used to construct a change of variables which 



will simplify the equations: we replace a and b by functions a(r) and 6(r), and set 

5^1 = 5<(a(r),6(r)), (57) 

6ij2 = 6^!^{a{r),b{r)) + ^2, 

5Z, = 5Z^{air),bir)) + Ze, 

5Z, = 5Z^{a{r),b{r)), 

5(P = 50"'(a(r),6(r)) + $, 

with 5K = 0. The first term on the right-hand sides is the diffeomorphism-induced 
perturbation (|39D for the diffeomorphism (^Pl), but with a and b now functions. 



Since the diffeomorphism satisfies the equations of motion for arbitrary constants a 
and b, the linearised equations will only involve derivatives of a(r) and b{r). The two 
first-order equations (p4| , p5D can then be solved for dra{r) and drb{r). Inserting these 



values into the other four second-order equations (p0|-p3|) and the dilaton equation 
(^) gives two equations which are trivially satisfied, and a coupled set of three 
second-order equations for "^2, and $. 

It is convenient to write the coupled equations so that each one only involves 
second derivatives of one of the functions. Then the equation which involves $" is 
(where ' again denotes 9^, and we assume that all the perturbations behave as e*'^*) 

^('^" + ^7^^' + ^^''^) + ^2 + 2^6) + + ^2^2 + Q\Z, = 0, (58) 

with the polynomial coefficients 

D = r^fC(8r^ + 5rr2 + Srfg + 2f2fe)(4r^ + 3rf2 + 3rfe + 2r2'f^6), (59) 
P2 = — 2r^fC(— 2r^f2 + Gr^fg + 8rf2fe + Sr^rQ + ^27^5), (60) 

Q\ = _r2(4rof2 + 36rof6 + 3f2 + 6f2f6 + 27f2) (61) 
— r(40ror2f6 + 2f2r6 + 30f2fg) — 12rorlfQ — Srlf^, 

Ql = ro{-2r'^f2 + Qr'^TQ + 8rf2f6 + Sflfe + f2fg), (62) 
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Ql = r2(8ror2 + 24rof6 + 9f2 + lOfafg + 9f2) (63) 
+r{2Arof2fe + Gflfe + 10f2^6) + ^^orlf^ - 2ror2fg. 
The equation involving is 

^(^2 + ^^'^2) + + PlZ', + + Ql^>, + QlZ, = 0, (64) 
where D is as before, and the other polynomial coefficients are 

P| = 64r^ + r^(-32ro + 120f2 + SSfg) (65) 
+r^(-76ror2 - UrofQ + 30fl + 172f2fe + 30fg) 
+r (— 15ror2 ~ 118rof2f6 — Ibrofg + AAf2'^6 + 52f2fg) 
+r(— 28ror2re — 36ror2fQ + Sf^rl) — ArQ-flrl, 

P| = -Sr^hKiSr^ + 16rre + Shh + 5f|), (66) 

Ql = 4f2(r^(-8ro - 6f2 - Qtq) + r{ArofQ - 7f2r6 + 3fQ) + 6rof2r6 + 2f2fg), (67) 

Ql = -2rof2 (Sr^ + I6rfe + + 5fl), (68) 

Ql = 4f2(r^(16ro + 18f2 + 2f6) + r(12rof6 + 21f2f6-fg)-3rof2f6 + 5rofg + 6f2fg). (69) 
The equation involving is 

^(^6 + ^^'^e) + P|^2 + ^1^6 + Q?*^ + ^2*2 + Ql26 = 0, (70) 

where D is as before, and the other polynomial coefficients are 

P2 = — 2r^if(6r^f2 — 2r^r6 + 8rr2re + f^re + 3f2fg), (71) 

p| = 64r^ + /(-32ro + 64r2 + 96f6) (72) 
+r^(-20rof2 - 52rofg + 30f^ + 76f2fg + 30f^) 
+r (— 15ror2 — 22ror2r6 — 15rorg + 28r2r6 + 20r2rg) 
+r{-12roflfQ - 4ror2fg + 8f2rg) - Aroflfl, 

Ql = r^(12rof2 + l2rofQ + 9f2 + 10f2f6 + 9fg) + r(8ror2^6 + lOflfe + 6f2fg) - Aroflfe, 

(73) 

<52 = ''o(6r^f2 - 2r^fg + 8rf2fg + fgfg + 3f2fg), (74) 

Ql = -r2(24rof2 + 8rof6 + 27f^ + 6f2f6 + 3f^) (75) 
-r(24rof2fg + 30f2fg + 2f2fg) + 2roflfQ - 6ror2^6 ~ ^fgrg). 

Leaving aside the decoupled mode 6ip3, which will be discussed in the next section 
(and which we will find leads to no instabilities), we have now reduced the perturba- 
tion problem to these three second-order equations. The background whose stability 
we are mainly interested in addressing is the shell branch solution, so we will also 
need to formulate appropriate matching conditions at the shell. The determination 
of the matching conditions and the numerical investigation of the existence of suit- 
able solutions of the equations for negative u;^ will be the subject of a forthcoming 
companion publication . 
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5 Stability of the free wave equation 

We will now discuss the stability of the solution against perturbation by just turning 
on Sip3. We make the ansatz 6ip3{t,r) = \l/3(r)e*'^*. Then ( ^9] ) implies 



Z2ZQ 



+ W^^3 



0. 



(76) 



If we consider the horizon branch, we need simply look for solutions of this equation 
regular on the horizon and at infinity. For the shell branch, ([76|) applies for r > r^, 
and 



Ho Ha 



L' 
T 



dl^^ + ( - + - 1 9,^3 



+ 



-CJ^\E'3 







(77) 



applies for r < re. Since the shell does not couple to 5ipz^ the appropriate boundary 
conditions at the shell are that \E'3 and dr^'i are continuous. 

We translate this into a standard one-dimensional bound state problem, by intro- 
ducing the tortoise coordinate 



Y WT:) Jr, I r < re 

Irl^df r>r,. 



This coordinate runs from — 00 at r = rg, to +cx3 at r = 00. We make a chan; 
variable^ ^ 

' -if) r > Tfi- 



Then (|77D becomes 
where for r > rg, 
W{r) = W> = 



ZoZp, 



tzl. 

4Z2 



(78) 

le of 
(79) 

(80) 



■4 [Z2 



1^' 




Z'2 




-- 
16 V 


Z2 


]- 


IK'' 




1 K 


^ rl{ 


5 



5 

16 



Z. 



1 Z2 Z'q 

8 Z2 Zq 



while for r < r^, 

W{r) = W< = 



L{r,) H2He 



IHl IHl 

4:H2 A He 



16 \H2 



16 \H, 



IH'^ 



8 H2 Hi 



6 



1 m\ IL^ 

4 \H2 hJ L r L 



12) 



'Note that 9r^3 being continuous is not equivalent to dr4' being continuous. 
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Plugging in the functions from (12), we have 



^> = T7r^^{[8(^2 + r6) + 16ro]r-3 + [3f2 + 3f2 + 30f2f6 + 20ro(f2 + f6)]r-^ 

+ [12f2f6(f2 + rg) + Qroih + fQ)'^]r~^ + [^flfl + 8rof2f6(f2r6)]^"'^ 
+4rof^f^r-^}. (83) 

The general form of iy< is complicated, but in the case M = 0, where we have simply 
an uncharged black hole inside the shell, 

= —'-^ . (84) 

Z2(re)Z6(re)L(re)r3 ^ ^ 

On the horizon branch, where r2 > 0, > everywhere, and there can be no 
instability associated with this mode. This is as we would expect; the horizon branch 
looks like a normal charged black hole solution, and the free wave equation has no 
non-constant solutions regular both on the horizon and at infinity. However, on the 
shell branch, there may be a region with < 0. (Since we take the horizon branch 
for the solution inside the shell, is always positive.) The leading term is always 
positive, as 

1 r-. ^ 1 



r2 + h + ro = -\l^rl + ± -^4r2 + rg > 0, (85) 
since \t2\ < tq. On the other hand, Wy is always negative near r = — f2- As r — > — f2, 

^ --4^ < 0. (86) 

If we considered just the pure repulson solution, this divergence would lead us to 
suspect the solution is unstable to a perturbation by dip^. Although one would need 
to consider the issue of boundary conditions at the singularity, diverges sufficiently 
quickly that there could be bound states supported away from r = — f2. The question, 
then, is whether the enhangon excises this instability, along with the various other 
undesirable features of the geometry. 

In figure we plot the potential for some representative values of the parameters. 
We see that there is a substantial region outside the shell where the potential is 
negative, and might suspect that this signals an instability. 

However, there is a general argument which says that there can never be an 
instability in this case |2^. First, we note that as ( [76l) is simply the free wave 



equation in this background, it always has dip'i = constant as a solution. In terms 
of the bound state problem (^OD , this translates into the statement that there is a 
zero energy [uj = 0) eigenmode ipo of (|5D|) which of the same sign and is bounded 
everywhere; we can take it to be always positive. This zero mode ipo does not vanish 
at the boundaries, so it is not a physical perturbation but it is still an acceptable 
mathematical solution of this equation. 

Now assume there is a discrete spectrum of bound states ipuj with negative energy. 
These are our hypothetical unstable modes with a;^ < 0. We can see from the form of 
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Figure 1: r^W plotted against t^/tq for (left) tq = lOrg, V = lOOOK, M = 0, (right) 
ro = re/lO, V = lOOOK, M = 0. 



(^) that they go to zero as r^, ±oo. This means that they are bounded solutions 
and physical perturbations of our problem. The standard 'node rule' for the number of 
nodes of the eigenf unctions of the discrete bound states says that in order of increasing 
energy, the nth eigenmode has n — 1 nodes (without including the boundary ones). 
Thus, the lowest negative mode ipuj^^^ must have no nodes in the sense of the above 
rule: we can take it to be everywhere positive. 

Both ipQ and ip^rr^ax solutions of the wave equation (pO|). By multiplying the 
equation for each mode by the other and taking the difference, and integrating over 
r*, we can obtain the equation 

/oo 
^CVr^axi^odn (87) 
-OO 

The left-hand side is the difference of the Wronskians calculated at the boundaries. 
Since the eigenmode ■j/'o approaches a positive constant at the boundaries = ±oo, 
while the eigenmode ipujmax goes to zero, the Wronskian vanishes at each boundary. 
Hence, the left-hand side is zero. On the other hand, since both ipui^ax ipo are sup- 
posed to be everywhere positive, the right-hand side cannot be zero unless Umax = 0. 
Thus, assuming the existence of eigenmodes ip^^ with tu^ < produces a contradic- 
tion. Hence there can be no such modes, implying that the geometry is stable to 
perturbation by Sip^. 

Acknowledgements 

We are grateful to James Gregory, Clifford Johnson and Matt Strassler for useful 
discussions and to Geoffrey Potvin for bringing to our notice a typo in (pl). AD is 
supported in part by EPSRC studentship 00800708 and by a studentship from the 
University of Durham. SFR is supported by an EPSRC Advanced Fellowship. 



18 



References 



G. T. Horowitz and R. C. Myers, "The value of singularities," Gen. Rel. Grav. 
27 (1995) 915-919, |gr-qc/9503062| . 



J. Polchinski and M. J. Strassler, "The string dual of a confining 
four-dimensional gauge theory," piep-th/0003136 . 



K. Pilch and N. P. Warner, "A^ = 1 supersymmetric renormalization group 
fiows from JIB supergravity," Adv. Theor. Math. Phys. 4 (2002) 627-677, 
[hep-th/0006066 . 



I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge 
theory: Duality cascades and xsb-resolution of naked singularities," JHEP 08 
(2000) 052, |hep-th/0007T9li 



J. M. Maldacena and C. Nunez, "Towards the large limit of pure A^ = 1 
super Yang Mills," Phys. Rev. Lett. 86 (2001) 588-591, jiep-th/OOOSOOl 



A. Buchel, A. W. Peet, and J. Polchinski, "Gauge dual and noncommutative 
extension of an A^ = 2 supergravity solution," Phys. Rev. D63 (2001) 044009, 
[hep-th/0008076 . 



N. Evans, C. V. Johnson, and M. Petrini, "The enhancon and N = 2 gauge 
theory /gravity RG fiows," JHEP 10 (2000) 022, |hep-th/000808li 



K. Behrndt, "About a class of exact string backgrounds," Nucl. Phys. B455 



(1995) 188-210, |Eii)-th/950610e . 

R. Kallosh and A. D. Linde, "Exact supersymmetric massive and massless 



white holes," Phys. Rev. D 52 (1995) 7137-7145, |hep-th/9507022 



C. V. Johnson, A. W. Peet, and J. Polchinski, "Gauge theory and the excision 
of repulson singularities," Phys. Rev. D 61 (2000) 086001, ^ep-th/991ll61 . 



D. Z. Freedman and J. A. Minahan, "Finite temperature effects in the 
supergravity dual of the A^ = 1* gauge theory," JHEP 01 (2001) 036, 
[hep-th/0007250 . 



A. Buchel, "Finite temperature resolution of the Klebanov-Tseytlin 
singularity," Nucl. Phys. B600 (2001) 219-234, |hep-th/00lH46| . 



A. Buchel, C. P. Herzog, I. R. Klebanov, L. A. Pando Zayas, and A. A. 
Tseytlin, "Non-extremal gravity duals for fractional D3-branes on the 



conifold," JHEP 04 (2001) 033, |hep-th/0 102105 



S. S. Gubser, C. P. Herzog, I. R. Klebanov, and A. A. Tseytlin, "Restoration of 
chiral symmetry: A supergravity perspective," JHEP 05 (2001) 028, 
liep-th/010"2172|. 



19 



C. p. Herzog and P. Ouyang, "Fractional Dl-branes at finite temperature," 
Nucl. Phys. B610 (2001) 97-116, ^-th/0104069| . 



S. S. Gubser, A. A. Tseytlin, and M. S. Volkov, "Non-Abelian 4-d black holes, 
wrapped 5-branes, and their dual descriptions," JHEP 09 (2001) 017, 
[hep-th/0108205 . 



C. V. Johnson, R. C. Myers, A. W. Peet, and S. F. Ross, "The enhangon and 
the consistency of excision," Phys. Rev. D 64 (2001) 106001, piep-th/0105077 



K. Maeda, T. Torii, M. Narita, and S. Yahikozawa, "The stability of the shell 
of d2-d6 branes in a = 2 supergravity solution," [hep-th/010706d| . 



M. Bertolini, T. Harmark, N. A. Obers, and A. Westerberg, "Non-extremal 



fractional branes," Nucl. Phys. B632 (2002) 257-282, |hep-th/0203064 



E. J. Weinberg, "Black holes with hair," |gr-qc/0106030[ 

M. Wijnholt and S. Zhukov, "Inside an enhancon: Monopoles and dual 
Yang-Mills theory," Nucl. Phys. B639 (2002) 343-369, |tiep-th/0 110109 



A. Dimitriadis, K. Maeda, S. F. Ross, and T. Torii, "Stability of the 
non-extremal enhangon II,". (To appear). 

S. S. Gubser and I. Mitra, "The evolution of unstable black holes in anti-de 



Sitter space," JHEP 08 (2001) 018, |hep-th/00lTT27 



R. Gregory and R. Laflamme, "Black strings and p-branes are unstable," Phys. 
Rev. Lett. 70 (1993) 2837-2840, |hep-th/930r05^ . 



C. F. E. Holzhey and F. Wilczek, "Black holes as elementary particles," Nucl. 



Phys. B380 (1992) 447-477, |hep-th/92020T4 . 



P. Kanti, N. E. Mavromatos, J. Rizos, K. Tamvakis, and E. Winstanley, 
"Dilatonic black holes in higher- curvature string gravity. II: Linear stability," 



Phys. Rev. D57 (1998) 6255-6264, |hep-th/9703T92 



20 



